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Abstract 

We study "draining bathtub" as an acoustic analogue of a three-dimensional rotating black hole. 
Rotating fluid near the sonic horizon necessarily gives rise to the superradiant modes, which are 
partially responsible for the thermodynamic quantities in this rotating vortex-like hole. Using the 
improved brick-wall method, we explicitly calculate the free energy of the system by treating the 
superradiance carefully and obtain the desirable entropy formula. 
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I. INTRODUCTION 



As suggested by Bekenstein, a black hole may have an intrinsic entropy proportional 
to the surface area at the event horizon^], and subsequently Hawking provided quantum 
field theoretic calculations for the Schwarzschild black hole|2[. Since then, there has been 
much attention to the statistical-mechanical origin of the entropy, especially for rotating 
black holes QDDDD . In the brick-wall method, quantum effects can be easily taken into 
account 0. Introducing a brick- wall cutoff makes it possible to remove the divergent term 
due to the infinite blue shift near the horizon The entropy from the brick- wall 

method consists of mainly two parts: the most dominant term compared to the logarithmic 
one gives the Bekenstein-Hawking entropy, and the other represents a typical infrared con- 
tribution at large distances. Although this original brick-wall method is useful for various 



models 



na na U4 



some difficulties may arise because it is 



assumed that there exists a thermal equilibrium between the black hole and the external 
field even in a large spatial region. Obviously, this method cannot be applied to a nonequi- 
librium system such as a system of non-stationary space-time with two horizons because 
the two horizons have different temperatures and the thermodynamical laws are also invalid 
there. Solving these problems, a thin-film method as an improved brick-wall method has 



been introduced 



19, 



. In the thin layer, local thermal equilibrium exists and the divergent 
term due to large distance does not appear any more. 

On the other hand, in Ref. many black hole issues have been treated as field the- 
oretical problems in fluid because this acoustic analogue was useful in that including its 
thermodynamics such as Hawking radiation and entropy might be tested hopefully in the 



laboratory. Moreover, a "draining 



black hole had been well defined 



bathtu 



22 



23 



}" referred to as an acoustic analogue of a rotating 



24 



. However, the conventional brick- 



wall cannot be applied in this model because it is impossible that the angular velocities 
of particles have a same constant value in whole region, while for the Banados-Teitelboim- 
Zanelli(BTZ) black hole|3 the method was able to be used to examine some results from 



superradiance effectively 



161 ] . Thus, in this paper, we would like to investigate the draining 



bathtub in terms of the thin-film method, which is helpful since the angular speed near 
the horizon can be approximately fixed to a constant. In Sec. [Hi the generic formulation 
of the free energy for a rotating black hole is given in the grand canonical ensemble, and 
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it will be shown precisely why the thin-film method should be used in our model. Then, 
in Sec. Hill the thermodynamic quantities are calculated by treating superradiant and non- 
superradiant (regular) modes more carefully. Finally, summary and discussions are given in 

Sec. us 



II. FORMULATION OF THE FREE ENERGY FOR A ROTATING BLACK HOLE 



In order to calculate the entropy of a given system in the original brick-wall method, 
we consider a quantum gas of scalar particles confined within a box near the horizon of a 
black hole and introduce a cut-off parameter as in Ref.j^]. The free scalar field is assumed 
to satisfy the Klein-Gordon equation given by (□ + m 2 )^ = with boundary conditions 

tfj(r H + h) = ?P(L) = 0, (1) 

where r# is the horizon, and rn + h and L represent the inner and outer walls of a "spherical" 
box, respectively, and h is an infinitesimal cutoff parameter. Suppose that this system is 
in thermal equilibrium at a temperature T = with an external reservoir. Then, in a 
stationary rotating axisymmetric black hole, a partition Z for an ideal bosonic gas in the 



grand canonical ensemble is given by 



16,Q 



lnZ = ]Tln]r[e-^-^f, (2) 

A fc=0 

where k = 0, 1,2, • • ■ is an occupation number, e\ and j\ denote the energy and angular 
momentum eigenvalues for a single-particle state A, respectively, and Q is the angular speed 
in equilibrium. The series in the partition function (J2J) has a finite value for e\ — Qj\ > 0, 
but it becomes divergent for e\ — Qj\ < 0, so it is ill-defined. In order to resolve such 
problem caused by the rotation of the geometry, we deal with the mode solutions of the 
Klein-Gordon field carefully, which will be of the form, ip(t,r,<p) ~ e -*wt+im<t>^ f or a Tes ^ 
observer at infinity(ROI) because there exist two Killing vector fields denoted by dt and d^. 

Note that the angular speed Q in Eq. (J2J) appears in the thermodynamic first law for a 
reservoir, i.e., dE = TdS + QdJ for a stationary rotating system. Besides, the angular speed 
of a particle for a ROI should be restricted because no particles can move faster than the 
speed of light. In fact, it takes a value between the maximum Q + and the minimum Q_ 
given by 



Ot(r) = flo(r) ± yj(d t ■ <V^ ■ d <t>? 1 d t • • fy, (3) 



where Qo(r) is the angular speed of Zero- Angular- Momentum-Observer fZAMO) jl^] . It is 
clear that both Q± converge to the constant value of Qh = ^o( r = t h) near the horizon, so 
the angular speed of every particles near the horizon can be always thought of as f2# . Since 
the dominant contribution to the physical quantities of the system, such as total entropy, 
is attributed to the quantum gas in the vicinity of the horizon, it is natural to assume that 
the system is in equilibrium with a uniform angular speed Q = Qh- 

Before formulating a generic free energy for a rotating black hole, the density function 
defined by g(cu, m) = dn{oj, m) / du is introduced in order to calculate the free energy strictly, 
where n(u,m) is the number of mode solutions whose frequencies, or energies, are all below 
ui for a given value of angular momentum m. Thus, g(u,m)duj represents the number of 
single-particle states whose energies lie between uj and u + du and whose angular momenta 
are m. Now, from the partition function (j2J), the free energy F is obtained as 

PF = - In Z = - E [ dug{u, m) In E [ e -# w - n * m) ] " . (4) 

It would be plausible to mention here that a ZAMO near the horizon(r m m) could measure 
only ingoing modes given by ip- m (x) ~ e -^t+«m<^ wm ] e a j^qj W ould see both ingoing and 
outgoing ones, where t = t, <f> = <f> — tint, Co = \uj — fi^m| > 0, and rh = sgn(co> — VLHm)m. 
Here, sgn(x) is 1 for x > and —1 for x < 0. The ingoing wave near the horizon consists 
of two parts; one is the so-called superradiant(SR) modes with lu — VLum < 0, and the 
other is the nonsuperradiant(NS) modes with u — VLh^l > 0. Then, e - lu)t + im <t> — ^t-im<t> 
for SR modes, and e - luJt + im ^ = e -^*+™^ f or jsjg mo des. Since only the ingoing modes are 
considered near the horizon, (e, j) has the value of (u,m) for single-particle states with the 
NS modes, while (e,j) becomes (—lu, —m) for the SR ones. Separating the SR modes from 
the NS ones, the free energy (UJ) should be replaced by F = F NS + F SR with 



PFns = E / m) In [l - e -^- Q " m A , (5) 

PFsr = E / d "9{", rn) In [l - e ^~^} . (6) 

Note that u is positive definite, and the density functions are given by g(u,m) = 
dn(u,m)/du for the NS modes and g(u,m) = —dn(u),m)/duj for the SR ones. Both 
Eqs. (0) and © can be obtained from (3F = — J2mf dug(u,m)hiJ2k ex P(~kpCj), where 
g(u, in) = dn/duj. 
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On the other hand, the angular speed of any particle cannot reach Qh over a critical radius 
in our model because of the restriction for Q from Eq. (jSJ), which will be explicitly shown in 
the following section so that global thermal equilibrium cannot be achieved. Instead, if we 
consider scalar particles confined within a thin layer near the horizon 19, 2(3], their angular 
speeds naturally take the same value of a constant Qh due to local thermal equilibrium, and 
the thin-layer method is useful to find the thermodynamic quantities in our model. 

Apparently, the degrees of freedom of a field within a thin layer near the horizon play a 
major role in the calculation of the entropy of a black hole; hence global thermal equilibrium 
is not necessary any more because particles are assumed to be distributed only in the narrow 
region. Since it is well known that Hawking radiation is derived from the vacuum fluctuation 
near the horizon, the Bekenstein-Hawking entropy should be associated with the field in this 
narrow region, where thermal equilibrium exists locally and statistical mechanics remains 
valid. This local thermal equilibrium is the main postulate of thin-film method, and the 
thermodynamic properties such as pressure and temperature near the horizon is assumed to 
vary slightly. The thickness of layer is supposed to be so small on a macroscopic scale that 
the physical quantities can be considered to be constant and that the narrow region can be 
locally in thermal equilibrium. Also, it is supposed to be very large on a microscopic scale 
to make sure that statistical mechanics remains valid. Specifically, in our model, the outer 
boundary L of spherical box in Eq. is replaced by r# + h + 5, where the parameter 5 is 
a positive physical small quantity related to the thickness of the layer. It means that S has 
the scale over Plank length, but the brick-wall cutoff h is very small compared to the Plank 
length. 



III. THERMODYNAMIC QUANTITIES 



We now set up an acoustic analogue of a rotating BTZ black hole in order to investigate 
its thermodynamics with superradiance taken into account. In the irrotational fluid, the 



propagation of sound waves is governed by an equation of motion 

1 



21] 



-9 



MV^sTd^) = o, 



(7) 
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where ip is the fluctuation of the velocity potential interpreted as a sonic wave function, and 
the metric is given by 



Po 

9^ = — 

c 



( 



-{c 2 -v 2 ) -vi 
J 

>lj 



with i,j = 1,2,3, (8) 



where c is the speed of sound wave, p an d v % Q are the mass density and the velocity of the 
mean flow, respectively, 5ij is the Kronecker delta, and Vq = S^v^Vq. Note that the velocity 
potential is linearized as \I/ = ipo + ip, and v = V"0o- 

We then consider a draining bathtub fluid flow described as a (2 + l)-dimensional flow 
with a sink at the origin. If the metric is stationary and axisymmetric, the equation of 
continuity, Stokes' theorem, and conservation of angular momentum yield that po is constant 
and ipo(i~, (ft) — —Aln(r/a) + B<p, where a, A, and B are arbitrary real positive constants 2^]. 
Then, the velocity of the mean flow is given by vq = —r(A/r) + <j){B/r). 

Now, let us consider the draining vortex case with A ^ 0. Dropping a position- 
independent prefactor from the metric (JHJ) , the acoustic line element for the draining bathtub 
is obtained as 

(A \ ^ / B \ ^ 
dr + — dt) + (rd(f) dt) , (9) 

where the radii of the horizon and the ergosphere are 



A yw+w 

r H = — , r e = , (10) 

c c 

respectively. However, the metric © makes it difficult to calculate thermodynamic quan- 
tities because of its (t,r)-component. Fortunately, this can be overcome by a coordinate 
transformation in the exterior region of A/c < r < oo. Using the transformation given 



by 



3, 27 



Ar AB 

dt -> dt + r 2 c 2 _ A 2 dr ' < ^~ > ^ + ^P) ' ir ' (U) 

the metric © can be rewritten as the conventional form, 

ds 2 = -N 2 dt 2 + N~ 2 dr 2 + r 2 (d(j) - tt dt) 2 (12) 

with 

N 2 (r) = l-^= r -^, no(r) = 4 = n H !f, (13) 

where we rescaled time coordinate by c for simplicity, and Qh = B/(cr 2 H ). Note that the 
metric (fT2*|) is similar to that of a rotating BTZ black hole, but two metrics have a little 



Q(r) The Angular Speed 




FIG. 1: The possible angular speed of a particle lies between the upper and lower curves, which 
denote the maximum and minimum angular speeds with respect to r, respectively. The shaded 
area represents the thin layer, which is located inside the ergoregion. 

difference: although setting J = 2B/c gives the same angular speed fio(r) of ZAMO, the 
lapse function N(r) has a different form from that of BTZ black hole, which is explicitly 
given by iV 2 = (r 2 — r 2 + ){r 2 — r 2 _)/{r 2 l 2 ). 

Now, the maximum and minimum angular speeds are obtained from Eqs. (J3J) and (|12|) as 

«±(r) = Oo(r) ± (14) 

As mentioned before, the angular speed of every particle is Qh in the vicinity of the horizon 
since Q± goes to f2# as r — > r#. Note that there exists a critical radius r c = \jr 2 H + fi^ 2 , 
where the maximum angular speed Q + is equal to Qh- This means that no particle can move 
along with the angular speed of Qh over r c as shown in Fig. ^ so the spherical box should 
be located inside the critical radius r c . Moreover, as discussed in the previous section, the 
radius of its outer boundary should be smaller than that of the ergosphere. Then, the radial 
part of the sonic wave satisfies 

+ r 2 iV 4 A; 2 (r; u, m)^ m (r) = 0, (15) 

where fc 2 (r;o;,m) = N~ 4 (uj — Q + m)(ou — f2_m). It can be easily shown that the function 
k(r; u, m) plays the role of the momentum eigenvalue in the WKB approximation. Therefore, 
in the thin layer with the range of rH + h<r<rn + h + 5, the discrete energy uj is related 



rN 



d_ 

dr 



rN 2 —i) u 
dr 



(r) 
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to the total number n(u,m) by 

7cn(uj, m) 



r H +h+S 



r H +h 



dr k(r; u, m), 



(16) 



where k(r;u,m) is set to be zero if k 2 (r;uj,m) becomes negative for given (u,m)8}. The 
contribution of k to the following calculations is dominant near the horizon because k is 
approximately N~ 2 and diverges as r goes to r#. This tells us that the thin-film method is 
valid. 

Then, we now evaluate the free energy of total system. The free energy for NS modes (0) 
is written as 

1 rr H +h+8 
7T Jr H +h 

k(r; u, m) 



/3F NS = W du— 

m Juj>Qj]in UUJ 



dr k(r; u, m) 



In 



-j3(io— Q[jm) 



f3 rr H 



r+h+5 



dr^2 / duj 



g/3(uj— Q^m) 



+ - 

7T Jr H +h 



dr J2 k(r- u, m) In [l - e -/?(«-n»»0 



(17) 



7T Jr H +h 
1 rrH+h+5 

I dr > /,•(/•: a.', m) lit 1 - i 

by using the integration by parts with respect to u. For the sake of convenience, the free 
energy Fns for NS modes is divided into two parts, which describe states with positive and 
negative angular momentum, i.e., F NS — ± NS -r j. ns 



Ft >0) +Fi™<°\ where 



-,(m>0) 



? (m<0) 
NS 



1 rr H +h+5 rao roo J (u) — Q + m) (oj — 0_77l) 

— / driV / dm / du;— ^ — = — r . 

7T Jr H +h JO Jn+m e ^-Q. H m) _ J 

1 pr H +h+S _ r 



(18) 



driV" 2 
X rm+h+s 



drN 



-2 



dm / duo 

oo JO 




JJoj — Q+m) (u — fl-m) 



II (3 Jr H +h 

Similarly, the free energy for SR modes (JOJ) becomes 



e (3(uj-n H m) _ X 

dmJQ + Q_m 2 ln(l - e ^ n * m ). 



(19) 



SR : 



r r H +h+s roo /-n_m «/ (u; — £7 + m) (a? — fi_m) 
drN / dm / du;— ^ — f; — ; 

7T Jr H +h Jo Jo e -P(u>-n H m) _ I 

I rr H +h+8 poo , 

— / drN~ 2 J dm^/n + n_m 2 \n(l-e- 



(20) 



7T/3 

Fortunately, the second terms of Eqs. (jTUjl and (|2Uj) are exactly cancelled, but large portion of 
tedious calculations is still required to be evaluated. After evaluating the above integrations 
with respect to u and m, the expression of total free energy is obtained as 



C(3) fTH+h+S 



4/3 3 Jr H +h 



dr- 



(q+ - n. 



N 2 (tt + - fi H )S(fi H -fi_)! 



(21) 



Then, substituting Eqs. (fT3|) and (JHJ) into Eq. (JUJ), the total free energy of our system is 
reduced to 

„ C(3H 



/3 3 



(22) 



2{h + 6) 

in the leading order. Note that there are no logarithmically divergent terms in the total free 
energy (|2*2l because those are remarkably cancelled as well as in the rotating BTZ black 
hole case[ly|. In addition, there exists no infrared divergence since the large distance is 
out of consideration, while infrared divergent terms remaining in total free energy were cut 
off from consideration in Refs.jsl, |^|. It seems appropriate to comment here that for the 
limiting case of non-rotating acoustic black hole, B = 0, there are no SR modes due to the 
fact that the angular speed of horizon vanishes, Qh = 0; in addition, there are no critical 
radius r c and no ergoregion due to r e = r#. Thus, one might think in this case having only 
NS modes the above result should be recast and give different value; however, the very same 
result is obtained so that Eq. (}22|) is remained valid in the limit of Qh — * 0. And since SR 
modes are distinguished from NS ones for ROI and not for ZAMO near the horizon, it is 
reasonable that the result from considering the superradiant modes in the rotating acoustic 
black hole is the same as the non-rotating one. 

On the other hand, the surface gravity is given by k 2 h = — ^VyV nXv 
where we used an appropriate Killing field near the horizon, x M = (d t + fi; 
then the temperature becomes 

1 



' r=rH r-l 1 /^' 

tf^ym and 



T, 



H 



Z7T 



2-Kr H 



Using the thermodynamic relation S = /3 2 dF/d/3\/3 = [3 H 
system is obtained from the free energy (f2*2~j) as 

3C(3) 



-3/3F| 



(23) 



the entropy of this 



S 



4vr 2 




2(h + 5) 



+ 0(Vh,Vh + 5) 



(24) 



and it can be rewritten in terms of the thin-wall cutoffs as 

3C(3) 6 



S 8tt 3 h(h + 5) 
where the cutoffs were defined as h = C H+h x / 

Jr H V 



+ 0(h,h + S), 



(25) 



r H +h+5 



g rr dr « y/2r H h and 5 = J^ +h 



2rn(h + 5) — y/2rnh in the leading order, and t = J ( 



2tt 



t=th 



g rr dr ps 
2ttth is the cir- 



cumference of the horizon. Note that h is called the brick-wall cutoff, and S represents 



9 



the thickness of the thin layer. Then, setting 5/[h(h + 5)] = 167r 3 /[3C(3)£ p ], or equiva- 



lently h = (5/2) [—1 + J 1 + 3C(3)£ p /(47r 3 5)], the entropy of sonic wave becomes finite and 
equivalent to the Bekenstein-Hawking entropy in the leading order, 

S = ^ = S B „, (26) 

where the three-dimensional Plank length is chosen as l v = hG/c 3 . It is plausible to make 
sure that the brick-wall cutoff becomes a universal value of h = [3C(3)/(167r 3 )]^ p in the 
leading order if 5 is larger than £ p , which is set to be one in the following calculations. 

Next, let us calculate the other thermodynamic quantities, such as the angular momentum 
of a matter particularly interpreted as a phonon in our model, 



dF 

•-'matter — 



= 3C(3y H Q H 5 
dn H ^ =fe 8tt 3 h(h + 8)' 1 ' 

where the partial derivative was evaluated from Eq. (j21j) . Note that substituting the expres- 
sion of cutoff h into Eq. ()27jl . the angular momentum is given by 

^matter = ^^H r H = ' 

The internal energy of the system in the frame of a ROI is explicitly calculated as 

E = F H + P^S + fWmatter = ^ + = 1 + (29) 

where Fh = F\p=p H . In the limit of non-rotating acoustic black holes of J ma tter — > 0, or 
equivalently B — > 0, it can be easily seen that the total energy E has the minimum value of 
4/3. 

Finally, it seems to be appropriate to comment on the perfect vortex case. If we take 
the limit of pure spinning acoustic black holes of A — > 0, the internal energy in Eq. (29) is 
undefined. Therefore, we must independently analyse this case whose spacetime represents 
a fluid with a non-radial flow. But, there does not exist the event horizon any more in this 
case. Therefore, it is meaningless to consider the analogy between gravitational and acoustic 
black holes. 



IV. DISCUSSION 



In this paper, we have studied the thermodynamics of a rotating acoustic black hole 
involving the superradiant modes for the draining vortex case (A ^ 0) as an acoustic analogue 

10 



of a black hole in three dimensions. In order to overcome some difficulties in applying the 
original brick-wall method to our model, the thin-film method has been introduced as an 
improved brick-wall one. And using this method we have obtained the thermodynamic 
quantities such as free energy, entropy, angular momentum, and internal energy of the thin- 
layer under thermal equilibrium with the black hole. The definition of thermodynamic 
black hole entropy was chosen in Eq. (}2l)j) as Sbh = 2£/£ p following that of BTZ black 
hole 22] to fix the brick-wall cutoff h, where the leading order of the entropy becomes S ~ 
Sbh{^ ~h/S) f° r a universal value of brick- wall cutoff. Recovering the physical dimension, 
the angular momentum (f2~%J) and the internal energy (|2TJ|) becomes J ma tter = 2c 2 B/G and 
E/c 2 = (4/3 + 2B 2 /A 2 )c 2 /G, respectively. 

As for the case of the limit of A — > 0, the internal energy (|29|) diverges. In fact, the 
metric (fT2*|) of A = seems to describe a pure rotation without horizons. However, this limit 
could not be taken since it has a naked singularity at r = 0. Therefore, we should consider a 
different form from the metric (fT2|) in order to deal with the pure rotation. Also, in the pure 
rotation, the brick-wall method could not be used to calculate thermodynamic quantities 
since the particles are distributed in whole region and it is impossible for the particles to fix 
the angular velocities to a special value. 

Finally, for the purpose of checking the stability of the system, the heat capacity can be 

n 

calculated as[3l| 

'dE s 



Cj= [dT, 



= 2S, (30) 

where we used the first law of thermodynamics, dE = TdS + QndJ, and the thermodynamic 
relation between the entropy and the free energy from Eq. ((2H)- Since the entropy is always 
positive, the heat capacity ()30|) is positive, and this means that the rotating acoustic black 
hole is thermodynamically stable. And also it can be easily shown that the curvature scalar 
of background geometry is positive everywhere as R = 2[r 2 H + (J/2) 2 ]/r 4 . 
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